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The J\f = 1 supersymmetric mKdV-B system is transformed to a system of coupled bosonic 
equations by using the bosonization approach. The bosonized supersymmetric mKdV-B 
(BSmKdV-B) equation includes the usual mKdV equation and a linear partial differential 
equation. The bosonization approach can thus effectively avoid difficulties caused by anti- 
commutative fermionic fields of the supersymmetric systems. The consistent tanh expansion 
(CTE) method is applied to the BSmKdV-B equation. A nonauto-Backlund (BT) theorem 
is obtained by using CTE method. The interaction solutions among solitons and other com¬ 
plicated waves including Painleve II waves and periodic cnoidal waves are given through a 
nonauto-BT theorem. The features of the soliton-cnoidal interaction solutions are investi¬ 
gated both in analytical and graphical ways by combining the mapping and deformation 
method. 
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I. INTRODUCTION 


The study of supersymmetric integrable systems has developed great importance in recent years 
00. The supersymmetrization of a number of integrable equations in which the bosonic equation 
is independent of the fermionic variable and the system is linear in fermionic field goes by the name 
B-supersymmetrization. It has received a lot of attention because of their connections with string 
theories from the point of view of matrix models |5|. The B-supersymmetric of the Korteweg-de 
Vries (KdV) 0], dispersionless two boson hierarchy [b], Sawada-Kotera [3], modified KdV and 
Camassa-Holme quations |8| have been constructed. In the meanwhile, the methodologies involved 
in the study of integrable systems have been expanded to the supersymmetric framework SB 
The soliton solutions of supersymmetric systems have been constructed via many methodologies. 
However, how to find exact interaction solutions among solitons and other kinds of complicated 
waves is an important problem for supersymmetric integrable models. In this paper, we will 
use the bosonization 14, 0 and consistent tanh expansion (CTE) 0 methods on the TV = 1 
supersymmetric mKdV-B (SmKdV-B) system. These interaction solutions among solitons and 
other types of solitary waves such as Painleve II waves and cnoidal waves are explicitly given. All 
results can be directly transformed to other B-supersymmetrization systems. 

The J\f = 1 SmKdV-B system reads [sj 


+ D 6 ® - 6 (D<$>) 2 D 2 <S> = 0, 


(1) 


where D = dg + 6d x is the covariant derivative. The commuting space variable x is extended to 
a doublet ( x,9 ), where 6 is a Grassmann variable. The Taylor expansion of the superfield with 
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respect to 9 is $(Q,x,t) = £(x,t) + 9u(x,t). The component version of (JX|) thus reads 

( 2 a) 
( 2 b) 


Ut 6 u U x T U xxx — 0, 
s/ 6u T Cxxx — 0. 


The bosonic equation (2a) is not depend on the fermionic variable, and equation (2b) is linear in 
the fermionic field. To avoid the difficulties in dealing with the anticommutative fermionic field of 
the supersymmetric equations, the component fields £ and u are expanded as the following form 
by introducing additional two fermionic parameters 14-16I. Il8|] 


^(x,t) = U 1 C 1 + U 2 C 2 , (3a) 

u(x,t) = U 0 + U 12 C 1 C 2 , (3b) 


where Cl an d C2 are two Grassmann parameters, the coefficients Uk = Uk(x,t)(k = 0,1,2) and 
U 12 = u\ 2 (x,t) are four usual real or complex functions with respect to the usual space-time 
variables. Substituting ([3]) into © yields the bosonized SmKdV-B (BSmKdV-B) system 


V'kjt 6Uq Uk t x T ^h.xxx 0, 

^ 12 , t 6 ( 1 X 12 ^ 0 )^ + Ui2,xxx — 0 . 


(4a) 

(4b) 


Equation (4a) is just the usual mKdV equation which has been widely investigated [19]. Equation 
(4b) is linear homogeneous in u\ 2 - These pure bosonic systems can be easily solved in principle. 
This is one of the advantages of the bosonization approach. 

The paper is organized as follows. In Section 2, a CTE approach is developed to the BSrnKdV- 
B equation. It is proved that the BSmKdV-B equation is CTE solvable system. A nonauto- 
BT theorem is obtained with the CTE method. In section 3, some novel exact solutions of the 
BSmKdV-B equation are derived through a nonauto-BT theorem. The last section is a simple 
summary and discussion. 


II. CTE METHOD FOR BSMKDV-B SYSTEM 


The CTE method is developed to find interaction solutions between solitons and any other types 
of solitary waves [l7j]. The method has been valid for lots of nonlinear integrable systems 2ol- 24]. 
According to the CTE method 2^25], the expansion solution has the form 


Uk = v = v 0 + v\ tanh(/), (5a) 

U 12 = w = wq + w\ tanh(/) + W 2 tanh(/) 2 . (5b) 


where vo, «i, u>o, wi, W 2 and / are functions of (x,t) and should be determined later. Substituting 
(5a) into (4a) gets 

6f x vi[y\ - / 2 ) tanh(/) 4 + 6(2 f x v\vQ + f£v 1:X + v 1 f x f xx ~ vjv 1>x ) tanh(/) 3 - (^f t + f xxx - 8/ 3 + 

6 f x vj - 6 f x vl + Sylvan + I2v 0 vi t x + ? ' (L ' 1 ^ x)x ^j tanh(/) 2 + + v hxxx - Q{vlvi) x - Qf x f xx vi - 

6/ 2 fi,x - 12 f x vfv 0 ) tanh(/) + v 0)t + v 0 , xxx + vi(f t + f xxx - 2/ 3 ) + 3(ui , x f x ) x - 6v%(v 0 , x + v x f x ) = 0. (6) 
By setting the coefficients of tanh(/) 4 and tanh(/) 3 in (0 to zero, one determines v\ and vq as 


Vl = ~fx , 


( 7 ) 
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By setting the coefficient of tanh(/) 2 in ([6]) to zero and using dT]), one obtains the equation for / 

ft = —fxxx + 2 f x + (8) 

^ Jx 

One can verify that the coefficients of tanh(Z) 1 and tanh(/)° in ([6]) are identically zero by using 
©. Substituting © and © into (4b) yields 

I2f x (flwi + f x w 2 , x ~ fxxW 2 ) tanh(/) 4 + 2(6 pw 0 - 9f x f xx wi + (2 f xxx - f t + 8/ 3 + 

~3f x W2,xx ) tanh(/) 3 - (3/ xx (2/ x + ^ - jf)w 2 + (ft - 5 f xxx + 10/ 3 - + 18f x f xx w 0 

^2 ,t ^2,xxx 4“ (24/j. 4 j/f)^ 2 ,z 3f xx wi jX + 3/ ' x wi iXX 4~ Qf x w o !X ) tanh(/) + (6 (f xxx 2/ 3 )ido 

+3f xx (2f x - fjr- + jt)w i 4- (2f t + 2/ xxx - 16/ 3 - ^f L )w 2 + wi,t + wi,xxz - (6/ x + ^£jr)wi, x 
+6(f x W2,x) x ~ 6f%wi,x) tanh (/) + 3f xx (2f x - ^ + jf)w 0 4- (ft + fxxx ~ 2 f x - 

+§f x f xx W 2 4- w 0 , t + w 0}XXX + 3(f x w hx ) x + 6 f x w 2:X - '' fx yi’ x = 0- (9) 


Similarly, by setting the coefficients of tanh(/) 4 and tanh(/) 3 in © to zero, one determines w\ 
and w o as 


W\ = 


f xx w 2 w 2tX 
P fx 


( ft f 

Wo = I- 


V6/ 3 3 /| 4/4 


+ 4% - iU + 


2/1 


3fxx^2,: 


( 10 ) 


We denote u /2 as n for simplicity. Setting the coefficient of tanh(/) 2 to zero, we reduce the equation 
for n 

_ rf2 _ 6/ xx T^XX H - 3fxxx^X 27f xx^X ^-2fXX fXXxTl !2 fln_ 

Tit i Tlxxx vj x Tlx p r /»3 v-L-v 

Jx *Jx Jx 

The coefficients of tanh(Z) 1 and tanh(/)° in © are identically zero by using dill) . While all the 
coefficients of powers tanh(/) of the BSmKdV-B system can be vanished by using appropriate 
Do, vi, wo, w i , w 2 and /, we call the expansion © is a CTE and the BSmKdV-B system is 
CTE solvable [17]. In summary, we have the following a nonauto-Backlund (BT) theorem for the 
BSmKdV-B system ©. 

Nonauto-BT theorem. The fields v and w 


ft 


v = —f x tanh(/) + —, 


, tfxxn n rx , , , „ s 4 n 


w = ntanh(/) 2 + (-g- - -f-) tanh(/) - -n + 

Jx Jx 3 *J: 


xx 

2 


(12a) 

9 fxx^x 4" 2 nfxxx 4~ Tift 3f xx n ^2b) 


3 ' 2 fl 6/ 3 2// 

are a solution of the BSmKdV-B system (f4j). while / and n satisfy © and (fill) . 


III. INTERACTION SOLUTIONS OF BSMKDV-B SYSTEM WITH A NONAUTO-BT 

THEOREM 

Some novel solutions of the BSmKdV-B system can be found by using the above nonauto-BT 
theorem. Here three examples are listed in the following. 

Example I. A quite trivial solution of © and dill) has the form 


/ = k\x + + li 


n = k 2 x + u 2 t + l 2 


(13) 
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where Aq, Aq, h and I 2 are arbitrary constants and oj\, ui 2 are determined by the relations 


uq = 2 kf, 


0J2 = Qk 2 ki. 


(14) 


The soliton solution of BSmKdV-B system reads in the following form by using the line solution 
m and the nonauto-BT theorem 


v = —ki tanh(Aq:r + 2 kft + l\), 

w = rotanh(Aqx + 2 k\t + l \) 2 — tanh(Aqx + 2 kft) — n. 

«i 


(15a) 

(15b) 


Though the soliton solution (11511 is a traveling wave in the (x,t) for the boson field v, it is not a 
traveling wave for other boson field w except for the case n being constants, i.e., k 2 = 0 . 

Example II. For the usual mKdV system, there exists a Painleve II reduction if one uses the 
scaling symmetry 26|. The scaling group invariant solution / of (JSJ) thus reads 


f = cln(t-b) + W(0, 


Z = 


x — a 


(16) 


(t-6) 3 

Substituting (fT6l) into (12a), a second order ordinary differential equation for the field W\ yields 


, 1 3WA 

= 2W 1 + + 


- C, 


W\ = W^, 


(17) 


which is the equivalent Painleve II reduction [17]], 0,6 and c are constants. The solution for other 
field n is obtained by solving 


m K£ = 


If + 6W?) M ( + WuiMj + WuMu 

27 W\ X M X + 12W hX Wi,xxM 12W^ X M 


(18) 


2 IF, 2 


Wf 


n = M( 0 , 


The interaction between solitons and Painleve II waves of the BSmKdV-B system can be obtained 
from the nonauto-BT theorem 


v = 


Wl - tanh(/) - 1 


VFl 5 


(t-b) 


2 {t-b)iWi 


(19a) 


w = M tanh(/) — 


2 M ? VFi - MW U 


wl 


tanh(/) - -M + 


9 - £M 
18 Wl 


^ cM - 9VFi i? M ? - 2MWi& 5MW^ 


6Wl 


4 Wf 


(19b) 


Example III. To find the interaction solutions between solitons and cnoidal periodic waves 
of (HI), the solution for the / field assume with one line solution k\x + u\t plus an undetermined 
traveling wave F(kx + ut) 


Substituting PUD into 


/ = k\x + uj\t + F(X), 
yields 

I2k \, 


X = kx + uit. 


Fi 2 x ~ 4Fi 4 - (Cik 3 + -^)Ff - (2>C\k\k 2 + 


12kf 2 lo 
~F r + 


)Fl 


( 20 ) 


( 21 ) 
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_ (3 Cl kkl + 4t? + U ' 1 + ^)*i - C!*f - Mfai + fcm) = 0 


k A 


k 5 


-^1 = 


with Ci is arbitrary constant. The solution of (12111 can be solved out in terms of Jacobi elliptic 
functions [271] - The solution expressed by (1201) is thus the explicit interaction solutions between 
one soliton and cnoidal periodic waves. The field n is obtained by solving the following equation 


Nxxx — 


6 k 


kF\ + k\ 


FixNxx + ( 


27 k 2 


3k 


+ 


UJ 

¥ 


N x +[ 


6k 2 


(kF 1 + k 1 ) 2 


Fi 


\2(kFi + k\) 
2k 


xx 


Fl ’X kF. + k/ 1 ^- 

k 2 

y F l,x) F l, X N - 0, 

n = N(X). 


( 22 ) 


As the well known exact Jacobi elliptic functions solutions of (l2lTl . we try to build the mapping 
and deformation relationship [151] between the solution for F\ and N by using (12211 . In order to get 
the mapping and deformation relationship, the variable transformation is introduced [151. Il6] 


Substituting the transformation 


N{X) = N\F 1 {X)). 

into (1221) and vanishing F\ x via (121 [1 . equation 


(23) 

becomes 


L F + §*(* + + 3M(Fl + - 3 MX - 0. ,24) 


where 


L = 4k 4 Ff + k 3 {C\k 4 + 8k x )F 2 + 2^0^^ + 2 kk\ + u)F 1 + C x k\k b + ku i + k x u, 

M = k 3 {C x k 4 - 4k\)Ff + 2k{Cik 5 h - 4 kkf + 2u)F 1 + C x h^k{ - 4 kk\ + 3ku x + k x uj, 

for simplicity. The solution for (12411 can be directly obtained by using Maple since (1241) is the 
linear ordinary differential equation. The mapping and deformation relation is thus constructed 
by solving Oil) 

N = N' = C 2 kF? + (C 2 k i + C 3 k)F x + C 3 k x , (25) 

with C 2 and C 3 are arbitrary constants. 

Substituting (l20l) . (l23l) and (l25l) into (fl2l) . the corresponding solution of the BSmKdV-B system 
reads 

v = —{kF x + k\) tanh(/) - , (26a) 

2 kF\ + k\ 

w = (kFi + k x )(C 2 Fi + C 3 ) tanh(/) 2 - C 2 kF 1:X tanh(/) + ^ 4C>2 ^^ 2 + F i,xx 

(26b) 

k 3 (3C 2 kF x + 2C 2 k x + C 3 k ) 2 (C 2 F X + C 3 )(8k 3 Ff + 24 k x k 2 F? + (24 kk\ - co)F 1 + 8 k\ - uj x ) 

4 {F x + fci)3 1>X 6[F\ + hY 2 

To show more clearly of this kind of solution, we give two special cases. 

Case I. One solution of m take as 


F\ = a-o + aisn^A - , m), 


(27) 
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FIG. 1: (a) Plot of one kink soliton on the cnoidal wave background expressed by (29a). (b) Plot of the first 
special soliton-cnoidal wave interaction solution by (29b). The parameters are k = 1, m = 0.3, a-2 = l,Ci = 
1, C'2 = 1 , C s = 1 . 


where ao, a\ and a 2 are constants, sn(a 2 Jf, m), dn(a 2 X, m) and cn(a 2 X,m) are the usual Ja¬ 
cobi elliptic functions with the modulus m. S , D and C stand for sn(o 2 X, m), dn(a 2 Jf, m) and 
cn(a 2 Jf, m) respectively for simplicity. Substituting (1271) into (l 2 TTl and vanishing all the coefficients 
of different powers of S, a nontrivial solution is obtained for constants 


oq — 


u = 


Ci A : 3 3a 2 

4 2~ 

k 3 a 2 {m 2 — 5) 


a 1 = 

uji = 


Q^vn 


= 


k(C\k 3 + 802 ) 


2 ’ 4 

k 3 a 2 (32a2 + 5CiA : 3 — C\k 3 m 2 ) 


(28) 


In this case, the interaction solution of of the BSmKdV-B system becomes by using 


v = —^ 02(1 + mS) tanh(/)--- —CD, (29a) 

2 2(1 + mS) 

w = -02^(1 + mS)(2C2Ci2mS + 4C3 — 6C202 — C1C2A: 3 ) tanh (/) 2 -f— C2C^2kmCI) tanh(/) 

8 2 


1 


+ -a 2 fc(CiC 2 A : 3 - 4C 3 )(1 + mS) + -a 2 2 C 2 kmS + 


a\C2k{m 3 S 3 — 2 m 2 + 3) 


4(1 + mS) 


(29b) 




The Hgure 1 plots the interaction solutions between solitons and cnoidal waves for the held v and 
w respectively. The parameters are k = 1, m = 0.3, a 2 = 1, k\ = 1, Ci = 1, C 2 = 1, C 3 = 1. 

Case II. Another special solution of (12111 reads 


F\ 


, , , _M 2 _ 

1 2 1 — nsn(&i X, m) 2 ’ 


(30) 


where b 1 , b 2 ■ n and m are constants, sn( 6 iX, m), dn( 6 iX, m) and cn( 6 iA, m) are the usual Jacobi 
elliptic functions. Here, S , D and C stand for sn( 6 iA, m), drv(biX,m) and cn( 6 iX, m) respectively 
for similarity. Substituting ( 13011 into (| 21 ll and vanishing all the coefficients of different powers of S 
yields constraints of constants 


h = 


C\k 3 -\/ n(l — n)(m 2 — n ) 
4 (m 2 n — 2m? + n) 


b2 = 


(1 — n)(m 2 — n ) 


n 


ki = 


Cik 4 (l — n)(m 2 — n) 
A(m?n — 2 m? + n) 


(31) 
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FIG. 2: (a) Plot of one kink soliton on the cnoidal wave background expressed by (31a). (b) Plot of a 
second special soliton-cnoidal wave interaction solution by (31b). The parameters are k = 1, m = 0.8, n = 
0.2, Cl =3,(72 = 1; C3 = 1. 


UJ = 


C\ k 9 (l — n){n — m 2 )(m 2 n — 3 m 2 + n) 


u 1 = 


Cfk 12 (m 2 n - m 2 — n 2 + n) 2 {m 2 — m 2 n — n) 


8 {m 2 n — 2m 2 + n) 2 ’ 32 (m 2 n — 2m 2 + n) 3 

In this case, a interaction solution for the BSmKdV-B system is obtained as 
C*iA; 4 (n — l)(m 2 — n) 


v = 


w = 


. ,/C nCi/cV(m 2 -n)(l-n) 

4(mW2m 2 + n)(nS 2 fl) ’ 

=/2 


4(n5 2 — l)(m 2 n — 2 m 2 + n) 

_/ CiC 3 fc 4 (n - l)(m 2 - n) nCfC 2 k 7 (n-l) 2 (m 2 -n) 2 - 2 

\4(nS 2 — l)(m 2 n — 2m 2 + n) 16(nS 2 — l) 2 (m 2 n — 2m 2 + n) 2 / 

+ CfC 2 k 7 n(m 2 - n)§(l - n)§ gg>£, tanll ^ + Cl<^3^(1 - n){m 2 - n) 


(32a) 

(32b) 


+ 


8 (nS 2 — 1 ) 2 (m 2 n — 2 m 2 + n) 2 4 (nS 2 — 1 )(m 2 n — 2 m 2 + n) 

C 2 C 2 k 7 n(n — 1 )(m 2 — n)(m 2 n 2 S 6 — 3 S A m 2 n + S 2 m 2 n + ^m 2 + 5 2 m 2 + S' 2 ??. — n) 


( / 


16(??S 2 — 1) 2 (?77. 2 ?? — 2m 2 + ??) 2 
C'f?n 2 A; 12 (l — n) 2 (m 2 — n) 2 t j (l — n)(m 2 — n) 


16(?n 2 n — 2 m 2 + n) 3 


n 


■E n (S, n, ??? 


where is the third type of incomplete elliptic integral. The figure 2 plots the interaction solutions 
between solitons and cnoidal waves expressed by ()32f) with the parameters k = 1, m = 0.8, n = 
0 . 2 , C\ = 3 , 6*2 = 1,63 = 1. There are some typical nonlinear waves such as interaction solutions 
between solitary waves and cnoidal periodic waves in the ocean [27F291]. Solutions (1291) and 
may be useful for studying these types of ocean waves 


28 


29(]. 


Remark. For a given solution v and w of BSmKdV-B system, the interaction solutions for 
component fields £ and u of SmKdV-B can be constructed by using ([3|). The bosonization approach 
can thus effectively avoid difficulties caused by anticommutative fermionic fields for supersymmetric 
nonlinear systems. 


IV. CONCLUSIONS 

In summary, the Af = 1 SmKdV-B system is mapped to a system of coupled bosonic equations 
by means of the bosonization approach. The BSmKdV-B system is just the usual mKdV equation 
together with a linear differential equation. Then, the CTE approach is applied to the BSmKdV-B 
equation. It is proved that the BSmKdV-B equation is CTE solvable system. A nonauto-BT 
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theorem is obtained with the CTE method. Various explicit solutions of the BSmKdV-B system 
such as soliton-Painleve II waves and soliton-cnoidal waves are obtained by using the nonauto-BT 
theorem. For the interaction soliton-cnoidal waves, two cases are given both in analytical and 
graphical ways via combining the mapping and deformation method. This kind of interaction 
solutions may be useful in real physical phenomena. 

For the bosonization approach, we can also introduce N > 2 fermionic parameters Ci (* = 
2,3, ...,1V) to study the B-supersymmetrization of nonlinear evolution system. On the other hand, 
lots of B-supersymmetrization systems are introduced from the usual action principle The 
interaction solutions among solitons and other complicated waves for these B-supersymmetrization 
systems are worth studying. 
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